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Final  Technical  Report  1994-1997 
AFOSR  F49620>96-l-0053 

**Compnter  Simulations  of  Growth  Phenomena  on  Solid  Surfaces** 

(a)  List  of  Objectives 

Our  work  provides  an  understanding  of  epitaxial  growth  at  the  atomic  level  and  develops  new 
phenomenological  theories  describing  epitaxy  on  large  space  and  time  scales.  There  are  practical 
reasons  for  seeking  such  an  understanding.  Solid  state  electronics  is  moving  relentlessly  toward 
smaller  and  smaller  devices  to  increase  both  storage  density  and  computing  speed  and  to  reduce 
power  consumption.  As  devices  get  smaller  and  the  surface-to-bulk  ratio  diminishes,  the  quality  of 
the  interface  becomes  very  important.  A  rough  interface  degrades  conductivity  and  luminescence 
and  increases  energy  loss.  Until  now,  progress  in  the  growth  of  small  and  complex  solid  state 
structures  has  relied  on  empirical  research.  It  is  hoped  that  a  better  understanding  of  crystal  growth 
will  allow  us  to  optimize  device  performance  and  perhaps  suggest  how  to  grow  new  structures, 
with  new  functions. 

(b)  Status  of  the  Research  Effort 

We  could  grow  a  perfect  solid  if  we  could  grow  a  perfect  layer  on  the  substrate  and  then  grow  the 
solid  one  layer  at  a  time.  This  is  why  a  lot  of  theoretical  and  experimental  effort  goes  into  studying 
sub-monolayer  growth.  Our  results  in  this  direction  arc  described  below. 

(1)  A  new  mechanism  for  surfactant  action.  Many  experimental  smdies  have  shown  that 
one  can  use  “surfactants**  to  improve  the  quality  of  epitaxially  grown  films.  These  findings  are 
especially  important  in  the  case  of  heteroepitaxy  where  lattice  mismatch  often  leads  to  rough  and 
fractured  films.  It  is  very  likely  that  surfactant  action  has  several  mechanisms  that  depend  on  the 

chemical  nature  of  the  system.  We  have  proposed  and  tested^  one  such  mechanism,  which  is 
operative  when  the  surfactant  atoms  repel  those  of  the  material  being  grown.  This  repulsion 
changes  die  character  of  the  atomic  motion  from  diffusion  to  a  slower  motion  that  we  call 
“anomalous  diffusion”.  This  in  turn  increases  the  number  of  islands  nucleated,  as  compared  to  the 
case  when  no  surfactant  is  present,  which  leads  to  layer-by-layer  growth.  We  established  this 
mechanism  by  analysis  and  computer  simulations,  and  have  shown  that  it  does  explain  the 
experimental  observations.  This  is  a  new  concept  in  cpitaXy,  likely  to  also  be  useful  in  the  kinetics 
of  surface  reactions. 

(2)  The  kinetics  of  dimers  in  the  early  stage  of  epitaxy  on  Si(  100).  In  the  past  several  years 
we  have  m^ped  the  kinetic  pathways  followed  by  the  first  Si  atoms  depo.sited  on  a  Si(lOO)  surface 
and  indicated  how  our  results  can  be  tested  by  low-tenqierature  STM  measurements.  Recent 
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developments  in  technology  have  made  such  measurement  possible.  They  confinned  our  early 

predictions  as  well  as  those  made  prior  to  the  experiments.^  We  have  predicted  that  the  Si  atoms 
deposited  on  the  surface  will  adsorb  on  top  of  the  dimer  rows,  will  move  very  rapidly  along  them, 
and  will  meet  and  make  dimers  and  longer  one*dimensional  clusters  on  top  of  the  rows.  At  higher 
temperatures  the  deposited  atoms  will  move  between  the  rows.  We  have  also  predicted  that  the 
dimers  on  top  of  the  rows  will  rotate.  These  predictions  were  rather  radical  at  the  time  they  were 
made.  The  diffusion  rate  given  by  us  was  much  larger  than  that  inferred  from  experiments  and  the 
clusters  we  predicted  were  oriented  perpendicular  to  the  ones  observed  at  high  temperature. 
Nevertheless  all  our  predictions  were  confirmed  by  detailed  direct  experim^ts. 

We  have  also  analyzed^  the  collective  properties  of  missing  dimers  on  Si(lOO)  and 
predicted  that  they  will  form  wavy  lines  on  the  surface.  These  have  been  detected  by  ZandvUet  in 
Holland  and  their  properties  agree  with  those  calculated  by  us. 

(3)  Island  shapes,  evaporation,  and  coarsening.  The  adsorbate  islands  formed  on  surfaces 
have  very  intciguit^  shapes.  The  equilibrium  shape  is  a  thermodynamic  quanti^  that  depends  on 
the  "surface”  tension  along  the  island’s  border.  In  principle,  by  analyzing  these  shapes  we  can 

calculate  this  tension.  We  have  developed  the  methodology^"5‘il  of  doing  so  in  two  cases:  when 
the  shape  is  determined  by  transport  of  material  along  the  border,  and  when  it  is  dominated  by 
evaporation/condensation  on  the  island.  We  have  shown  in  both  cases  how  one  could  use  the 
measured  shape  to  determine  the  edge  tension.  The  evaporation  of  atoms  from  the  islands  leads  to 
"coarsening”:  a  change  in  island  density  and  size  distribution.  In  catalysis  this  is  a  highly 
obnoxious  event  which  degrades  the  performance  of  the  catalyst  and  needs  to  be  avoided.  In 
microelectronics  it  is  a  desirable  effect  since  it  can  be  used  to  make  island  sizes  more  uniform.  We 
have  developed  new  methods  to  simulate  coarsening  and  in  the  process  have  studied  the 

ev^>oration  from  islands  in  detail.^*!  l  We  discovered  that  the  evaporation  rate  has  strange 
properties  and  explained  their  origin.  We  have  also  shown  that  island  coalescence  follows  a 
scaling  law  *3  which  we  have  tested  by  simulations. 

(4)  Growth  of  heterogeneous  systems.  A  system  is  heterogeneous  if  the  lattice  of  the  top 

layer  is  mismatched  with  that  of  the  substrate.  A  good  example  is  Au  on  a  Au(lOO)  surface, 
Instead  of  making  a  square  lattice  (as  would  be  expected  on  a  (100)  substrate),  the  top  Au  layer  has 
a  nearly  hexagonal  packing  and  its  density  is  20%  higher  than  that  of  the  bulk.  The  islands  grown 
on  this  stuface  are  rectangular  and  show  a  "size  quantization”:  they  tend  to  be  six.  or  twelve,  or 
eighteen  atoms  wide.  We  have  used  an  effective  medium  theory  to  calculate  the  activation  barriers 
for  atom  motion  on  and  around  the  islands,  to  discover  the  reasons  for  the  quantization  of  size  and 
the  preference  for  long  rectangles.  The  bond  length  between  the  gold  atoms  on  the  top  layer  is 
shorter  ttian  in  the  bulk.  The  atoms  tend  to  get  close  to  each  other  in  the  horizontal  direction  and  to 
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A  Directed  Algorithm  for  Finding  Potential  Energy  Functions,  Materials  Research  Seminar, 
Institute  for  Theoretical  Physics,  l|niversity  of  California,  Santa  Barbara,  May  27, 1997 
Searching  for  Potentials  and  Bootstrapping  Time  Scales,  UCSB  Institute  of  Theoretical  Physics 
Workshop  on  Interatomic  Potentials  and  Linking  of  Scales,  June  5-7, 1997 
Activated  Dynamics  on  Surfaces,  Euroconference  on  Technical  Advances  in  Particle-Based 

Computational  Material  Sciences  “Computer  Simulation  of  Rare  Events  and  the  Dynamics  of 
Classical  and  Quantum  Condensed  Phase  Systems’*,  Lerici,  Italy,  July  7-1 8, 1997 
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buckle  vertically,  so  that  along  one  direction  six  atoms  in  the  island  fit  on  top  of  five  atoms  on  the 
substrate.  Because  of  this  buckling  the  mobility  of  an  atom  across  an  island  is  low,  and  the 
mobility  along  an  island  is  high.  This  leads  to  growth  of  long  islands.  Moreover,  the  diffusion 
constant  of  an  atom  along  the  long  edge  of  a  rectangular  Island  depends  on  the  size  of  the  island 
and  is  fastest  if  the  island  is  six,  twelve,  or  eighteen  atoms  wide.  Once  an  island  reaches  one  of 
these  widths  the  atoms  reaching  its  border  move  rapidly  along  the  long  edge  of  the  rectangle,  and 
reach  the  short  edge  and  get  stuck  there.  This  mechanism  favors  both  the  long  rectangle  shape  and 
size  quantization.  We  are  convinced  that  diis  mechanism  ($ize*dependent  mobility  across  die  island 
or  along  its  edges)  is  essential  in  understanding  growth  in  heterogeneous  systems. 

{Sy  Thermodynamic  and  kinetic  studies  of  the  mean  island  shape  and  its  surface  tension. 
The  thermodynamics  and  the  kinetics  of  the  small  islands  grown  on  a  solid  surface  are  fascinating 
because  of  the  smallness  of  the  system.  If  we  think  of  such  an  island  as  a  large  molecule,  which  it 
is,  then  STM  pictures  of  evolving  islands  provides  us  with  information  on  the  isomerization  of  a 
molecule  in  contact  with  an  energy  reservoir.  Thermodynamics  has  been  developed  to  deal  with 
the  mean  properties  of  a  huge  number  of  molecules.  What  is  then  the  meaning  of  thermodynamics 
when  applied  to  measurements  that  study  one  island?  The  next  question  is  how  to  implement 
thenTOdynamics  to  study  a  system  in  which  all  measurements  are  kinetic  and  we  cannot  measure 
the  two  essential  quantities,  work  and  heat  We  have  addressed  this  question  in  two  lectures  at  the 

NATO  school*®  in  Rhodes.  We  proposed  that  the  thermodynamic  theory  is  meaningful  only  if  it 
deals  with  averages  over  an  ensemble  of  islands.  We  showed  how  such  an  ensemble  can  be 
constructed  by  performing  kinetic  measurements  or  simulations.  This  has  been  applied,  in 
conjunction  with  kinetic  Monte  Carlo  simulations,  to  deteinune  the  surface  tension  of  the  island’s 
edge  from  its  equilibrium  with  a  gas  of  atoms  adsorbed  on  the  substrate,  and  also  from  the 
evolution  of  its  shape.  The  results  of  these  two  independent  theories  are  consistent  with  each 
other. 

(6)  A  thermodynamic  theory  of  evaporation  rate.  In  work  sponsored  by  this  grant,  we 
made  a  surprising  discovery:  the  rate  of  evaporation  of  atoms  from  an  island  is  proportional  to 

N®-^,  where  N  is  the  number  of  atoms  in  the  island.  The  exponent  was  independent  of  all  other 
variables  in  the  system,  temperature,  size,  and  kinetic  model.  Everyone  in  the  field  expected  the 
exponent  to  be  1/2.  To  understand  this  size  dependence  we  developed  a  theory  of  the  evaporation 

9 

rate  which  predicts  that  the  evaporation  rate  is  proportional  to  N*^  exp[  b/N^^].  This  is  even 
more  surprising  and  seems  to  contradict  our  earlier  results.  However,  both  functional  forms  (i.e. 
NO.36  cxpl  b/N^^J )  fit  the  data  equally  well!  The  constant  b  in  the  second  equation 

depends,  in  a  known  way  ,  on  the  edge  tension  of  the  island.  By  performing  simulations  of 


F49620-96-1-0053 


3 


nHK-ix-X3:?a 


ucoo  urriwc  ur  ni.coc.nrs.<^n 


0:7*J  CUXX 


r  • 


ev^qporation  rates  we  extracted  b  and  the  surface  tension.  The  latter  agrees  very  well  with  the 
surface  tension  extracted  from  simulations  of  island-vapor  equilibrium'^,  proving  that  our  theory  is 
sound. 

(7)  Aspects  of  the  coarsening  process.  If  we  grow  an  ensemble  of  islands  and  they  do  not 
have  the  desired  size  or  position  distribution,  we  have  one  more  chance;  by  heating  the  system  we 
can  cau.se  the  atoms  of  the  islands  to  evaporate,  move  on  the  substrate,  and  join  other  islands. 
Heating  also  causes  the  islands  to  move  along  the  surface,  meet  each  other,  and  merge.  This 
process,  known  as  coarsening,  will  change  the  island-size  distribution  and  the  way  the  islands  are 
arranged  in  space.  We  have  studied,  by  Kinetic  Monte  Carlo  simulations,  the  mobility  of  the 
islands  and  the  change  in  island-size  distribution  caused  by  coarsening.  It  was  proposed  in  the 
literature,  by  many  audiors,  that  the  diffusion  coefficient  of  a  migrating  island  has  a  power-law 
dependence  on  the  size,  with  a  universal  exponent  (i.e.  independent  of  material  or  working 
conditions).  Our  simulations  show  that  this  is  not  the  case.  Recent  experiments  in  Wendelkcn’s 

group  have  confirmed  our  findings.  We  have  also  developed  a  phenomenological  theory  of  the 
coarsening  process  taking  place  during  deposition,  as  the  islands  grow  in  size,  touch,  and  merge. 
We  found  a  new  scaling  law  for  island-size  distribution  and  confirmed  it  by  simulations. 

(8)  The  theory  of  rate  constants.  Kinetic  simulations  require  the  computation  of  a  large 
number  of  rate  constants.  Activated  complex  theoiy  is  a  favorite  method  for  determining  the  rates. 
Such  calculations  have  two  difHcult  steps:  finding  the  transition  state  and  convincing  oneself  that 
the  transition  state  approximation  is  conect  We  have  proposed  a  new  procedure  for  finding 
transition  states.  This  used  a  genetic  programming  method  to  “teach”  the  ab  initio  program  the 

configurations  for  which  to  calculate  energy  to  find  the  transition  state  most  efficiently.'’  We  have 
also  studied  a  variational  method  for  finding  the  best  dividing  surface  for  the  transition  state,  in 
complex  systems  mat  have  multiple  recrossings.'^ 
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18.  Island  Migration  Caused  by  the  Motion  of  the  Atoms  at  the  Border;  Is  There  a  Scaiing  Law?, 
A.  Bogicevic,  S.  Liu,  J.  Jacobsen,  B.  Lundqvist,  and  H.  Metiu,  Phys.  Rev.BIS  (Rapid 
Communications),  to  appear 
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19.  Thermodynamics  and  Kinetics  of  Two-dimensional  Islands  on  Surfaces:  A  Simulation 
Study,  Lone  M0llnitz,  Master’s  Thesis  (1997),  Niels  Bohr  Institute,  University  of 
Copenhagen  (supervised  by  H,  Metiu  and  J.  Norskov) 

20.  Migration  of  Hydrogen  on  a  Solid  Surface:  The  Physics  of  the  Process  and  the 
Methodology,  H.  Metiu,  Proc.  Fermi  School  on  the  Dynamics  of  Rare  Events  (Lerici,  July 
1997) 

2 1 .  Some  Constraints  Involving  the  Statistical  Properties  of  Trajectories  Run  in  the  Monte  Carlo 
Computation  of  a  Rate  Constant  and  Their  Use  in  Improving  and  Testing  the  Quality  of 
Sampling,  D.  E.  Makarov  and  H.  Metiu,  J.  Chem.  Phys„  to  appear 

(d)  Professional  personnel  supported  by  the  grant 

Mr.  Alexander  Bogicevic  (visiting  from  Chalmers  University),  Dr.  Lutz  Bdnig,  Mr.  John  Detch 
(undergraduate),  Dr.  Shudun  Liu,  Dr.  Dmitrii  Makarov,  Dr.  Gregory  Mills,  Ms.  Lone  M0llnitz 
(visiting  from  the  University  of  Copenhagen),  Dr.  Hongxiao  Shao,  Dr.  Seokmin  Shin,  Dr.  Paul 
Weakliem 


(e)  Interactions 

Presentations  at  seminars,  meetings,  conferences 
Seminars 

Crystal  Growth,  Island  Shapes  and  Nucleation,  Technical  University  of  [.amannft,  December  1995 
Crystal  Growth,  Island  Shapes  and  Nucleation,  University  of  Geneva,  December  1995 
Electrons  in  Zeolites,  Centre  d’Etudes  Nucleaires  de  Grenoble,  December  1995 
The  Properties  of  Electrons  Solvated  in  Zeolites,  Tel  Aviv  University,  January  1996 
Adsorbate  Self-Aggregation:  Kinetics  and  Thermodyruxmics,  Joint  Chemical  Physics-Theoretical 
Physics  Seminar,  Fritz  Haber  Institute,  Berlin,  March  7, 1 996 
The  Properties  of  Electrons  Solvated  in  Zeolites,  Freie  Universitat,  Berlin,  March  8, 1996 
Adsorbate  Self-Aggregation:  Kinetics  and  Thermodynamics,  Physics  Department,  Aarhus 
University,  April  11, 1996 

Adsorbate  Self-Aggregation:  Kinetics  and  Thermodynamics,  Physics  Department,  Technical 
University  of  Denmark,  Lyngby.  April  1 6, 1996 
The  Properties  of  Electrons  Solvated  in  Zeolites,  University  of  Paris,  Sorbonne,  May  6,  1996 
Lectures  at  Professional  Meetings 

The  Electronic  and  Magnetic  properties  of  Akali-Doped  Zeolites,  International  Chemical  Congress 
of  Pacific  Basin  Societies,  Honolulu,  December  1995  (presented  by  Nick  Blake) 
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An  Investigation  of  the  Electronic  and  'Optical  Properties  of  Sodium  Clusters  in  Sodalite, 

International  Discussion  Meeting  on  Time-Dependent  Quantum  Mechanics  of  Many-Electron 
Systems.  Bangalore.  January  9-12, 1996 

Adsorbate  Self-Aggregation,  Symposium  on  Reaction  Dynamics  at  Surfaces,  Orsted  Institute, 
Copenhagen,  March  1, 1996 

Electrons  in  Zeolites,  Plenary  Lecture  at  the  Annual  Meeting  of  the  German  Zeolite  Society  ,  Berlin, 
March  5, 1996 

The  Electronic  and  Magnetic  Properties  of  Alkali-Metal-Doped  Sodalites,  ACS  National  Meeting, 
New  Orleans,  March  1996  (presented  by  Nick  Blake) 

Kinetic  Simulations  of  Island  Formation  During  Epitaxial  Deposition,  AFOSR  Contractor’s 
Meeting,  Colorado  Springs,  June  2, 1996 

An  Investigation  of  the  Electrordc  and  Optical  Properties  of  Sodium  Clusters  in  Sodalite,  ACS 
National  Meeting,  Orlando,  August  25, 1996 

Two  lectures  at  the  NATO  ASI  Workshop  on  “Surface  Diffusion:  Atomistic  and  Collective 
Processes”,  Rhodes,  August  27-Septcmber  6, 1996 

(a)  Beating  the  Time  Scale:  Simulations  of  Femtosecond  Events  All  the  Way  to  Minutes  of 
Real-Time  Runs 

(b)  Diffusion  Controlled  Aggregation  of  Adsorbed  Atoms:  Island  Densities.  Shapes, 
Evaporation,  and  Coarsening 

Lectures  on  Surface  Processes:  (i)  Beating  the  Time  Scale:  Simulations  of  Femtosecond  Events  All 
the  Way  to  Minutes  cf  Real-Time  Runs,  (ii)  Diffusion-Controlled  Aggregation  of  Adsorbed 
Atoms:  Island  Densities,  Shapes,  Evaporation,  and  Coarsening ,  NATO  ASI  Workshop  on 
Surface  Diffusion:  Atomistic  and  Collective  Processes,  Rhodes,  August  27-September  6, 
1996 

Lectures  on  Zeolites:  (i)  Zeolites:  Structure  and  Uses,  (ii)  Zeolites:  Interaction  Potentials.  (Hi) 
Isolated  Electrons  Solvated  in  Zeolites,  (iv)  The  Properties  of  Fully  Doped  Zeolites: 
Semiconductors  or  Metals?,  Australian  Summer  School  on  Condensed  Matter  Physics, 
Canberra,  January  10-31, 1997 

Simulations  of  Epitaxial  Growth,  Seminar  Lecture,  Qiemistry  Department,  Hong  Kong  University 
of  Science  and  Technology,  Febmary  3, 1997 

Initial  Stages  of  Epitaxial  Growth:  Modeling  and  Simulations,  Syn^)Osium  on  Epitaxial  Growth, 
Materials  Research  Society  Meeting.  San  Francisco.  March  31-April  4, 1997 
Thermodynamics  and  Kinetics  of  Small  Adsorbate  Islands  Studied  by  KMC  Simulations,  ACS 
National  Meeting,  San  Francisco,  April  13-17, 1997 
Simulations  of  Epitaxial  Growth,  Materials  Science  Seminar,  Applied  Physics  Department, 
California  Institute  of  Technology,  May  14, 1997 
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